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Double perovskites like Sr2FeMo06 are materials with half- metallic ground states and ferrimag- 
netic Tc's well above room temperature. This paper is the second of our comprehensive theory for 
half metallic double perovskites. Here we derive an effective Hamiltonian for the Fe core spins by 
"integrating out" the itinerant Mo electrons. We next validate the classical spin Hamiltonian by 
comparing its results with those of the full quantum treatment presented in the companion paper 
"Double Exchange Mechanism for Half-metallic Double Perovskites". We then use it to compute 
magnetic properties as a function of temperature and disorder. We discuss the effect of excess Mo, 
excess Fe, and anti-site disorder on the magnetization and Tc. We conclude with a proposal to 
increase Tc without sacrificing carrier polarization. 



Strong electron correlations and the interplay among 
charge, spin and lattice degrees of freedom lead to a 
wide range of spectacular phenomena in transition metal 
oxides^ such as high Tc superconductivity, colossal mag- 
netoresistance and large thermopower. 

Half metals with fully spin polarized ground states pro- 
vide another example of such unique and spectacular phe- 
nomena. Among the known examples, double perovskites 
(DPs) are of particular interest due to their high ferro- 
magnetic Tc's along with the possibility of integrating 
different functionalities with oxide electronics^. One of 
the best-studied half-metallic DP is SrsFeMoOe (SFMO) 
with Tc=420 K, well above room temperature^^. DPs 
have the form A2BB'06 , which is derived from the sim- 
ple ABO3 perovskite structure with a three-dimensional 
(3D) checkerboard ordering of B and B' ions. DPs have a 
range of fascinating properties from spin liquids to mul- 
tiferroics, as well as from metals to multi-band Mott 
insulators ^^i^i^ — 

This article is the second part of our comprehensive 
theory for half metallic double perovskites. We begin 
by summarizing the first paper (hereafter referred to as 
paper I) titled "Double Exchange Mechanism for Half- 
Metallic Double Perovskites"^ where we discuss the full 
quantum Hamiltonian describing core spins on Fe cou- 
pled to conduction electrons through a generalized double 
exchange mechanism. In paper I, we calculated the mag- 
netic and electronic properties as a function of tempera- 
ture using exact diagonalization of the "fast" electronic 
degrees coupled to "slow" core spins configurations gener- 
ated by classical Monte Carlo simulations (ED+MC). By 
retaining the electronic degrees of freedom, we obtained 
information about the temperature dependent density of 
states and the destruction of the fully polarized half- 
metallic ground state through thermal fluctuations. One 
of our central results is that the conduction electron po- 
larization at the chemical potential is directly propor- 
tional to the core-spin magnetization. This finding is 
significant because it indicates that if one can derive an 
effective Hamiltonian for the core spins, it would be pos- 
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FIG. 1: (a) Schematic showing energy levels at transition 
metal sites in two unit cells (formula units) of SFMO. The Fe 
sites have localized 5* = 5/2 core spins, treated as classical 
vectors with orientation {0,(j)). The parameters t,t' and A 
of the Hamiltonian, Eq. ([!]), governing the dynamics of the 
itinerant electrons in t2g orbitals, are also shown, (b) Nearest 
neighbor (NN) and next nearest neighbor (NNN) configura- 
tion of two unit cells of DPs. 



sible to deduce the electronic polarization, a quantity of 
central importance for spin injection and spin transport, 
but one that is difficult to measure directly. A second 
problem the effective Hamiltonian will help ameliorate is 
the possibility of simulations on large system sizes com- 
pared to severe size limitations faced by ED+MC meth- 
ods. 

With this motivation, here we derive the effective 
Hamiltonian (Heff) for the Fe core spins by integrating 
out the itinerant Mo electrons. This is a non-trivial gen- 
eralization of the previous Anderson-Hasegawa analysis 
for manganites^ now applied to double perovskites. We 
validate Heff by comparing the temperature-dependent 
magnetization M(T) curve with that obtained from the 
ED+MC method for the full Hamiltonian. We find that 
the modified double exchange form of Heff indeed cap- 
tures the magnetic properties of the full Hamiltonian at 
all temperatures; retaining only Heisenberg interactions 
works only at intermediate temperatures. It is remark- 
able that Heff also captures possible metallic antiferro- 
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magnetic instabilities with increasing carrier concentra- 
tion. Regarding the effects of disorder, while both excess 
Fe and Mo decrease the saturation magnetization and Tc, 
anti-site disorder in which Fe and Mo exchange places, 
behaves differently; although magnetization drops, Tc is 
not affected. This finding forms the basis of our proposal 
to increase Tc without sacrificing conduction electron po- 
larization. Our results suggest that by putting excess Fe 
and compensating the loss of carriers with La doping can 
indeed lead to a dramatic increase in Tc- 

We start by briefly describing the full quantum Hamil- 
tonian. We then solve the problem of two unit cells and 
derive the effective exchange Hamiltonian between two Fe 
core spins. We then generalize this form to the infinite 
lattice. 

For large Hund's coupling Jh, Fe^+ in the 3d^ config- 
uration saturates the "up" manifold and forms a large 
spin S=5/2 that we treat classically with a local axes of 
quantization along S^. Mo^+ (4d^) contributes to con- 
duction in t2g orbitals. Due to the symmetry of t2g or- 
bitals, dai3 orbitals can only delocalize in a/3 planesi^ 
{a/3 = xy^yz^xz). For all the Mo sites j, we choose the 
same (global) axis of quantization. The generalized dou- 
ble exchange Hamiltonian^ that describes the core 
spins interacting with conduction electrons is 

H = -t ^ {dadl^Cja + h.c.) 

I 



-t' E ^lcr-^^T.4id^l (1) 

UJ'),cr i 

where dio- (Qct) are fermion operators on the Fe (Mo) sites 
with spin a. 

The orientation (^^, ^i) of the classical spins affects 
the Mo-Fe hopping via e^^ = — sin(^i/2) exp(i^i/2) and 
= cos(6>i/2)exp(-z0i/2). 

I. EXACT SOLUTION OF TWO SITE 
PROBLEM 

Next we solve the Hamiltonian, Eq. ([T]) exactly ana- 
lytically for two unit cells, shown schematically in Fig. [TJ 
This is a generalization of the Anderson and Hasegawa 
analysis for manganites^ applied to double perovskites. 

In a single unit cell, there are three states derived from 
the Fe^ and Mo^^^ t2g orbitals. We label the unit cells as 
i and j, and without loss of generality choose a coordi- 
nate system such that one of the core spins is aligned 
with the z axis, and the other core spin Sj lies in the x-z 
plane (Fig. HJa)). This particular choice of coordinates 
simplifies the calculation as it gauges away the depen- 
dence. Thus, = sm{6i/2) and = cos(^i/2) , where 
is the relative angle between and Sj. The two unit 
cell Hamiltonian is given by 
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in the basis of {Fe^, Moi^, Moi^, Fej^, Moj^, Moj^}. Here 
7=1 for nearest neighbor (NN) and for next nearest 
neighbor (NNN) configurations (See Fig. mb)). By con- 
verting the 6x6 matrix for H in a block diagonal form, 
it can be solved analytically. 

The eigenvalues are only a function of the angle be- 
tween the core spins Si and Sj and describe the effective 
magnetic exchange Hamiltonians. The lowest eigenvalue 
for NN configuration, describing one electron in two unit 
cells corresponding to n=0.5, is 

= - x/(A/2)2 + 2t2(l + cos(0/2)) (3) 

or equivalently 

if™ = + 2tHl + S, • S,)/2 ) (4) 

where S is the unit spin vector. We obtain a very inter- 
esting modified functional form with a double square root 



structure that is different from conventional Heisenberg 
or previously studied Anderson-Hasegawa models^. Note 
that the interaction is ferromagnetic and we calculate the 
spin stiffness JpM = d'^E/dO'^ by expanding the energy 
close to i9 = 0, where E{e) ^ E{0) + {l/2){d'^E/de'^)0'^ + 
0((9'^). We find, 

( -t fort » |A| 
\ -tVA, for t « |A| 

showing that the kinetic energy of the conduction elec- 
trons sets the scale of the ferromagnetic exchange. 

For two electrons in two unit cells, which corresponds 
to n=l, the effective Hamiltonian is obtained by adding 
up the lowest two eigenvalues. For the NN configuration, 
effective Hamiltonian is antiferromagnetic given by 

H^i = -v /(A/2)2 + 2^2(1 + cos(^/2)) 

-V(A/2)2 + 2^2(1 -cos(^/2)). (6) 
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II. MAGNETIC PROPERTIES OF THE 
EFFECTIVE HAMILTONIAN 

In this section we use the ferromagnetic Heff to calcu- 
late the magnetic properties of SFMO. This is justified 
because SFMO is deep in the ferromagnetic phase and 
phase competition is not relevant as discussed in paper I 
|8|- The electronic density of SFMO is n=0.33 and the 
lowest density achievable in the two unit cell calculation 
is n=0.5, we will nevertheless use the ferromagnetic Heff 
form of the NN and NNN exchange interaction generalize 
it to the infinite lattice 



FIG. 2: Energy as a function of 6 for ferromagnetic Heff (one 
electron in two unit cells) and antiferromagnetic Heff (two 
electrons in two unit cells) . Effective Hamiltonian gives hints 
for filling dependent magnetic phase transition. We include 
FM Heisenberg Hamiltonian (Hneis) for comparison. Note 
that FM Heff is quadratic for a broader range of compared 

to HHeis. 
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FIG. 3: (a) Spin wave spectrum of full Hamiltonian and the 
Heff, (b) M(T) comparison between full Hamiltonian, Heff, 
and Heisenberg Hamiltonian. All simulations are done with 
an 8 X 8 system due to the high computational cost of the exact 
diagonalization and Monte Carlo calculations. 



Upon increasing the electron density (n=0.5 1), we 
find that the effective magnetic coupling changes from 
ferromagnetic to antiferromagnetic which is rather un- 
conventional. The filling driven FM-AFM transition has 
also been discussed by others^^. Metallic antiferromag- 
netism with large local moments at a commensurate wave 
vector is rare in nature. The exchange stiffness is given 

by, 



Jaf 



t 

tVlA3| 



for t > |A| 
for t < lAI 



(7) 



with the scale for antiferromagnetism also set by the ki- 
netic energy of the conduction electrons. Even at a two 
unit cell level, Hgff provides a hint for this transition and 
illuminates the mechanism, though only discrete fillings 
are accessible at this level. 



^eff = -Jl ^1 • ^j) -J2YF2 (S, • S,) (8) 

(hj) {{hj)) 



The function Fi{x) = 8\/2 + + 2x is obtained by 
setting A=0 in Eq. (jH for the NN configuration and 
a similar procedure for the NNN configuration yields 
F2{x) = (5 + V^)\/6 + 2V3 + 2x. By expanding Heff 
for small 0, we obtain the spin wave spectrum in terms 
of the exchange constants Ji and J2, whose values can 
now be fitted by matching with the spectrum of the full 
quantum Hamiltonian, as shown in Fig. [3l^a). The agree- 
ment over the entire spectral range rather than just at 
small energies is indeed remarkable. The specific pref- 
actors (8 for Fi and (5 + Vb) for F2) are included for 
convenience to match the spin wave stiffness with that of 
the Heisenberg model with nearest Ji and next neighbor 
J2 interactions. 

Validation of the magnetic H^f^: We have performed 
ED+MC calculations for the full quantum Hamiltonian 
along with classical Monte Carlo simulations for Heff and 
for the Heisenberg Hamiltonian. In Fig. [3fb) we present 
a comparison of the temperature-dependent magnetiza- 
tion M(T) calculated for each of these three Hamilto- 
nians on an 8x8 system. It is encouraging to observe 
that M(T) calculated from Heff agrees remarkably well 
with that obtained for the full Hamiltonian at all tem- 
peratures^ whereas the Heisenberg Hamiltonian fails at 
intermediate temperatures T2^Tc/2. This result gives us 
confidence in the validity of Heff . 

Note that the full quantum model where classical 
spins are coupled to conduction electrons has low-lying 
fermionic excitations. If one thinks in a functional in- 
tegral language, integrating out the fermions would give 
rise to various extra exchange terms like longer range in- 
teractions and four or more spin exchanges. The fact 
that we can reproduce M(T) at all temperatures shows 
that the effect of such terms are negligible and we cap- 
ture the most important magnetic exchange interactions 
within our model. 

The agreement between M(T) for Heff and the full 
quantum Hamiltonian also indicates that both Ji and 
J2 are temperature independent. Although it is not clear 
a priori why this is the case, the fact that Tc is much 
less than the bandwidth provides a reasonable justifica- 



4 



4 

3 

1 






(a) 


■ oL=16 


\ T 


- n 12 




■ o8 

■ 1 ■ 





T<T^ 

c 




^(b) 










V 


T — 1 ^ 
J^— ID 


A 


L=12 




k o 


L=8 








'mil 1 







200 400 
T(K) 



eL 



50 



FIG. 4: (a) Magnetization as a function of temperature, 
M(T), of Heff by classical Monte Carlo calculations for in- 
creasing 3D system sizes: 8*^, 12^ and 16^. (b) Estimating 
the thermodynamic Tc using finite size scaling. M(T) for dif- 
ferent system sizes collapses to a universal function close to 
Tc with universal critical exponents. We used 3D 0(3) uni- 
versality class exponents and e = \T — Tc\/Tc is the reduced 
temperature. As a result, we found Tc=0.14t for SFMO. 
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FIG. 5: Types of disorder: (a) excess Fe, (b) excess Mo, (c) 
anti-site disorder. Black, blue and red lines represent FM 
bonds, the broken FM bonds and the superexchange between 
Fe sites. 



tion for the temperature independence of the exchange 
constants up to temperatures of order Tc. 
Phase Transition and Determination of Tc : The primary 
advantage of the classical Hamiltonian Heff is our ability 
to simulate much larger system sizes compared to those 
using ED+MC methods. We have performed the first 
3D finite temperature simulations of magnetic properties 
using classical Monte Carlo on up to 16^ unit cells on an 
FCC lattice, as shown in Fig.HJ^a)). We have determined 
Tc using the finite size scaling of M(T). According to 
the finite size scaling hypothesis, M(T) for a system of 
size is described by a function of the form M{T^ L) = 
L-^/^J'(eLi/^) where J^{x) is a universal function and 
e = |T - Tc|/Tc. The critical exponents p = 0.36 and 
v = 0.70 are known for the 3D 0(3) universality class. 
Using Tc as a fitting parameter, we plot M(e)L^/^ against 
eL^/^ for L = 8, 12 and 16. For the true thermodynamic 
Tc all curves, of different system sizes, collapse onto a 
single curve, as shown in Fig. IH^b) providing an estimate 
of Tc = 0.14t for SFMO. 

Low temperature spin wave contribution to M(T): Stan- 
dard ferromagnetic spin waves produce a T^/^ reduction 
of the magnetization, also known as the Bloch T^/^ law^^«. 
However, in Fig. [3fc), M(T) is linear at low T and this 
linear behavior in fact persists up to a relatively large 
fraction of Tc. We explain this difference, between the 
Bloch Law and the calculated linear behavior, as aris- 
ing from the difference between classical and quantum 
magnons. The classical Hamiltonian is equivalent to tak- 
ing the S ^ oo limit of the quantum Hamiltonian but 
keeping Tc ~ JS'^ constant. The T^/^ law is restricted 
to a temperature scale ksTo ^ Tc/S^ the magnon band- 
width or equivalent ly to Tq/Tc ~ 1/S. Therefore the 
range of temperatures to observe the Bloch law is com- 
pletely quenched in classical calculations and highly sup- 
pressed in the experiment due to the large S=5/2 on Fe. 

In order to understand the origin of the linear temper- 
ature dependence of the magnetization, we consider the 



role of spin waves to get 



M(T) = Mo 



l^tB.Z. 
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(9) 



where the integral is over the first Brillouin zone. For 
small the dispersion for magnons Wq ^ q^. As 5 ^ 
oo the exponential can be expanded at all temperatures: 
^IBJSwq _ g ^ 1 + ^Z^^ fQj. constant Tc. Upon 



using this expansion and evaluating the integral gives 
M(T) - Mo(l - 7T) where 7 = 0(1). Thus classical 
spin waves indeed provide a natural explanation of the 
linear T dependence of the magnetization at low T. 

The reason for the robustness of the the linear M(T) 
dependence up to relatively high temperatures is the 
peculiar double square root form of Heff (See Fig. [2j). 
Compared to the Heisenberg Hamiltonian, Heff is har- 
monic {E ~ JO'^) for a larger domain of 9. There- 
fore magnon-magnon scattering which is mainly due to 
the non-harmonic part of the Hamiltonian is highly sup- 
pressed and that explains why the spin wave regime and 
correspondingly the linear T behavior of M(T) survives 
up to relatively high T. Similar M(T) has been observed 
in experiments both on single crystals^^ and on thin 
films^^. 



III. DISORDER 

In SFMO, there are three common types of disor- 
der: excess Fe, excess Mo and anti-site disorder (See 
Fig. [5]). By using Heff, we perform large scale calcula- 
tions of the temperature-dependent magnetic properties 
on systems up to 16^ to investigate the effects of disor- 
der. Finite size effects close to Tc are highly suppressed 
with increasing system size as shown in Fig. HJ^a). We 
start the discussion with the general chemical formula 
Sr2Fei+2yMoi_2y06 with y greater (smaller) than zero cor- 
responding to excess Fe (Mo), followed by anti-site disor- 
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FIG. 6: Effects of Fe k Mo disorder for SrsFei+^Moi-^Oe 
using Heff and comparing it with experiments, (a) Fe rich 
(y>0) M(T), (b) Mo rich (y<0) M(T), (c) Tc as a function of 
y, (d) Saturation magnetization, M(0), with y compared with 



experiments 



der. We conclude with a proposal to increase Tc without 
sacrificing conduction electron polarization. 
Excess Fe: For y>0, as seen Fig. [5ja) Fe replaces Mo 
sites which has two main effects: First, it reduces the 
total conduction electron density that weakens the dou- 
ble exchange mechanism. Secondly, when two Fe sites 
are close to each other, the strong antiferromagnetic su- 
perexchange locks the spins. We estimate the strength 
of this superexchange S{S -\- 1) Jaf ^ 34 meV based on 
Tn=750K for a similar compound LaFeOs with S=5/2 
spins on Fe. The excess Fe spin with the down orienta- 
tion on the Mo site couples antiferromagnetically to the 
four neighboring up spins creating a local puddle that 
enhances ferromagnetism in its neighborhood. Capitaliz- 
ing on this enhanced ferromagnetism will form the basis 
of our proposal to enhance Tc- 

In Fig. ini^a), shows that the saturation magnetization 
M(0) drops with increasing amount of excess Fe, largely 
because of its antiferromagnetic coupling to the neigh- 
boring Fe sites (see Fig. [6l^d)). For small values of y, Tc 
does not change significantly, then drops rapidly (Fig. 
[6]^d)) beyond 0.05. The initial insensitivity can be 
attributed to the changes in Tc from the two effects, re- 
duction of conduction electrons and formation of local 
ferromagnetic puddles, balancing each other. The be- 
havior of both M(0) and Tc as a function of y are in 
good agreement experiments^^. 

Excess Mo: Excess Mo (y<0) leads to a dilution of 
the ferromagnetic bonds (see Fig. [5j) as well as an in- 
crease in conduction electron density. The detrimen- 
tal effects of dilution and broken ferromagnetic bonds 
on the magnetization as a function of T shown in Fig. 
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FIG. 7: (a) Anti-site disorder results for Tc and sat- 
uration magnetization, M(0) (both normalized with re- 
spect to their disorder-free values) compared with experi- 
ments, (b) Proposal to increase Tc by La and Fe doping, 
Laa;Sr2a;Fei+yMoiy06. Tc(y) for compensated (x=3y) and 
uncompensated (x=0).The uncompensated Tc(y) is compared 
with experiment^i^. 



ini^b) is reflected directly in the linear decrease of sat- 
uration magnetization M(0) and Tc as a function of y 
(see Fig. [6fd)). Once again these results are in good 
agreement with experiments^^. The behavior of M(0) in 
off-stoichiometric SFMO is also in agreement with DFT 
calculations^^. 

Anti-site disorder: An example of anti-site disorder (AS) 
in which Fe and Mo sites replace each others is shown in 
Fig. Efc). This is the most prevalent type of disorder in 
SFMO. It can be thought of as a combination of excess 
Fe and Mo disorder while keeping the carrier density con- 
stant. We quantify AS disorder using 5 the fraction of 
Fe atoms that are on the Mo sublattice; ^=0.5 is a fully 
disordered system. Fig. [Tl^a) shows that M(0) drops lin- 
early with a slope of (1-2^), primarily due to the flipping 
a Fe spin from the parallel to the antiparallel direction 
on a Mo site per AS pair, as shown in Fig. Efc). Tc ap- 
pears to be insensitive to AS disorder, primarily because 
two effects balance each other. While the broken EM 
bonds in the Mo rich regions weakens EM, the puddles 
of Fe rich regions has the opposite effect. Although Fe 
sites are coupled antiferromagnetically in these puddles, 
it locally creates stronger ferromagnetic domains. We 
believe that these two effects balance each other and Tc 
does not change significantly with anti-site disorder with 
very good agreement with experiments^^. 
Proposal to increase Tc : We conclude with a proposal to 
increase Tc without sacrificing conduction electron polar- 
ization. We propose adding excess Fe, that locally creates 
strong ferromagnetic puddles, and simultaneously adding 
extra La to compensate the loss of carriers. Our results 
are shown in Fig. [7Jb) and suggest that with adequate 
amount of La doping, Tc can be increased by about lOOK. 

The general formula for both La and Fe doping is 
Laa;Sr2-xFei+2yMoi_2y06. Assuming that the Fe valency 
remains fixed at +3, and only Mo valency changes from 
+5 to +5+7^ with doping, the charge balance dictates 
that r]={2y - x)/(l -y) or correspondingly, carrier con- 
centration changes from n=l/3 to n = (1 + x - 3y)/3. 
This implies that y=3x is the necessary amount of La 
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doping to keep n constant and the result of imposing this 
constraint is shown in Fig. El^b). Next we argue that our 
approach for enhancing Tc is better than only La doping. 
It is known that La substitution of x = 1 gives rise to a 
15% increase of Tc^^. However, this is accompanied by 
a huge increase in the extent of anti-site disorder The 
effect of La doping x=l, changes the Mo valence from 
+5 to +4 (using 77 = —x). The reduced electrostatic 
attraction between the Mo and the surrounding oxygen 
octahedra leads to an expansion of the MoOe octahe- 
dra. As the volume of the MoOe octahedra approaches 
that of FeOe, the B-B' ordering becomes fragile^ and 
the increased anti-site disorder reduces the polarization 
significantly^. In contrast, our proposal suggests a 25% 
increase in Tc is obtained for y = 0.25 and x = 0.75, with 
an average Mo valence of only +4.66 which is unlikely to 
give rise to large amounts of anti-site disorder. 

Finally, we have checked that the proposed system with 
excess Fe and La compensation is indeed fully polarized 
at T=0 using ED+MC. The increase in Tc by about lOOK 
is extremely encouraging as that would increase the room 
temperature polarization significantly. 



magnetism in the double perovskite half metals. The ef- 
fective magnetic Hamiltonian Heff with the double square 
root form, obtained after integrating out the itinerant 
electrons, is very different from standard Heisenberg or 
double exchange Hamiltonians and agrees remarkably 
well when compared with the full quantum Hamiltonian. 
Heff is found to retain the harmonic 0^ form in the cant- 
ing between neighboring spins up to a larger range of 
0. As a result classical spin waves provide a good de- 
scription of the temperature dependent M(T), with sup- 
pressed magnon-magnon scattering. We have performed 
large scale simulations of Heff with different types of dis- 
order. From our insights on the dependence of the sat- 
uration magnetization and Tc on disorder, we propose a 
mechanism to substantially increase Tc by balancing ex- 
cess Fe doping and compensating the loss of carriers with 
La doping. 
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